{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import sys\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "sys.path.append('../../')\n",
    "sys.path.append('../')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "PATH = \"/home/zionxie/Documents/completely/\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "from bin.utils import walk_dir,read_mat"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "路径 /home/zionxie/Documents/completely/ 扫描完成！\n"
     ]
    }
   ],
   "source": [
    "xp,yp = walk_dir(PATH)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "def draw_line(path):\n",
    "    label = read_mat(path,'y')\n",
    "    plt.plot(list(range(885)),label)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "def draw_hist(path):\n",
    "    label = read_mat(path,'y')\n",
    "    plt.hist(label,bins=10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAD4CAYAAADrRI2NAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAAPXklEQVR4nO3df6zdd13H8eeLlYHjVwe9LLPtvFUKOqeG5WaOkCAygmUj6xLJsgWkYGMjTkQgQoHEGQ3JFhSEBJHC5orBwZzoGgfiUkcWDZ3eMRj7wY869qO1oxdh88ciMHz7x/mSXEu7e+75nnPPbT/PR3Jzvt/P99c7n9y+7qefc77fk6pCktSOJ0y7AEnSyjL4JakxBr8kNcbgl6TGGPyS1Jg10y4AYN26dTU7OzvtMiTpuHLbbbd9s6pmlnvcqgj+2dlZ5ufnp12GJB1Xktw/ynFO9UhSYwx+SWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1Jglgz/J1UkOJ7nzKNvekqSSrOvWk+T9SfYnuSPJ2ZMoWpI0umFG/NcAW45sTLIReBnwwKLmlwObu58dwAf7lyhJGqcl79ytqluSzB5l03uBtwI3LGrbCny0Bt/usi/J2iSnV9WhsVQrNWR2541Tue59V1wwletq5Yw0x59kK3Cwqr54xKb1wIOL1g90bUc7x44k80nmFxYWRilDkjSCZQd/klOAdwC/2+fCVbWrquaqam5mZtnPGJIkjWiUh7T9BLAJ+GISgA3A55OcAxwENi7ad0PXJklaJZY94q+qL1XVs6tqtqpmGUznnF1VDwF7gNd0n+45F3jE+X1JWl2G+TjntcDngOclOZBk++Ps/ingXmA/8GHgN8ZSpSRpbIb5VM+lS2yfXbRcwGX9y5IkTYp37kpSYwx+SWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktSYUb56UZ3ZnTdO5br3XXHBVK4r6cTgiF+SGmPwS1JjDH5JaozBL0mNWTL4k1yd5HCSOxe1vTvJl5PckeSvk6xdtO3tSfYn+UqSX5pQ3ZKkEQ0z4r8G2HJE203AWVX1s8BXgbcDJDkTuAT46e6YP0ly0tiqlST1tmTwV9UtwLeOaPv7qnqsW90HbOiWtwIfr6rvVNXXgf3AOWOsV5LU0zjm+H8V+HS3vB54cNG2A12bJGmV6BX8Sd4JPAZ8bIRjdySZTzK/sLDQpwxJ0jKMHPxJXgu8AnhVVVXXfBDYuGi3DV3bD6mqXVU1V1VzMzMzo5YhSVqmkYI/yRbgrcCFVfXook17gEuSPCnJJmAz8M/9y5QkjcuSz+pJci3wYmBdkgPA5Qw+xfMk4KYkAPuq6ter6q4k1wF3M5gCuqyqvj+p4iVJy7dk8FfVpUdpvupx9n8X8K4+RUmSJsc7dyWpMQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1BiDX5Ias2TwJ7k6yeEkdy5qe2aSm5J8rXs9tWtPkvcn2Z/kjiRnT7J4SdLyDTPivwbYckTbTmBvVW0G9nbrAC8HNnc/O4APjqdMSdK4LBn8VXUL8K0jmrcCu7vl3cBFi9o/WgP7gLVJTh9TrZKkMRh1jv+0qjrULT8EnNYtrwceXLTfga7thyTZkWQ+yfzCwsKIZUiSlqv3m7tVVUCNcNyuqpqrqrmZmZm+ZUiShjRq8H/jB1M43evhrv0gsHHRfhu6NknSKjFq8O8BtnXL24AbFrW/pvt0z7nAI4umhCRJq8CapXZIci3wYmBdkgPA5cAVwHVJtgP3Axd3u38KOB/YDzwKvG4CNUs6Qc3uvHEq173vigumct1pWTL4q+rSY2w67yj7FnBZ36IkSZPjnbuS1JglR/xS66Y1/SBNiiN+SWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BLUmOO+xu4vLlGkpbHEb8kNcbgl6TGHPdTPVpZPjZXOv4Z/JL+H983O/E51SNJjTH4JakxBr8kNcbgl6TGGPyS1JhewZ/kTUnuSnJnkmuTPDnJpiS3Jtmf5BNJTh5XsZKk/kYO/iTrgd8C5qrqLOAk4BLgSuC9VfUc4NvA9nEUKkkaj75TPWuAH0myBjgFOAS8BLi+274buKjnNSRJYzRy8FfVQeAPgQcYBP4jwG3Aw1X1WLfbAWD90Y5PsiPJfJL5hYWFUcuQJC1Tn6meU4GtwCbgR4GnAFuGPb6qdlXVXFXNzczMjFqGJGmZ+kz1vBT4elUtVNX3gE8CLwTWdlM/ABuAgz1rlCSNUZ/gfwA4N8kpSQKcB9wN3Ay8sttnG3BDvxIlSePUZ47/VgZv4n4e+FJ3rl3A24A3J9kPPAu4agx1SpLGpNfTOavqcuDyI5rvBc7pc15J0uR4564kNcbgl6TGGPyS1BiDX5IaY/BLUmP8zl0dF/weWGl8HPFLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUmF7Bn2RtkuuTfDnJPUlekOSZSW5K8rXu9dRxFStJ6q/viP99wN9V1U8CPwfcA+wE9lbVZmBvty5JWiVGDv4kzwBeBFwFUFXfraqHga3A7m633cBF/UqUJI1TnxH/JmAB+LMktyf5SJKnAKdV1aFun4eA0452cJIdSeaTzC8sLPQoQ5K0HH2Cfw1wNvDBqno+8N8cMa1TVQXU0Q6uql1VNVdVczMzMz3KkCQtR5/gPwAcqKpbu/XrGfwh+EaS0wG618P9SpQkjdPIwV9VDwEPJnle13QecDewB9jWtW0DbuhVoSRprNb0PP4NwMeSnAzcC7yOwR+T65JsB+4HLu55DUnSGPUK/qr6AjB3lE3n9TmvJGlyvHNXkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8kNcbgl6TG9L2BS1Mwu/PGaZcg6TjmiF+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfkhpj8EtSYwx+SWqMwS9Jjekd/ElOSnJ7kr/t1jcluTXJ/iSfSHJy/zIlSeMyjhH/G4F7Fq1fCby3qp4DfBvYPoZrSJLGpFfwJ9kAXAB8pFsP8BLg+m6X3cBFfa4hSRqvvs/j/2PgrcDTuvVnAQ9X1WPd+gFg/dEOTLID2AFwxhln9CxDkkY3ze+4uO+KC1b8miOP+JO8AjhcVbeNcnxV7aqquaqam5mZGbUMSdIy9RnxvxC4MMn5wJOBpwPvA9YmWdON+jcAB/uXKUkal5FH/FX19qraUFWzwCXAP1TVq4CbgVd2u20DbuhdpSRpbCbxOf63AW9Osp/BnP9VE7iGJGlEY/my9ar6LPDZbvle4JxxnFeSNH7euStJjTH4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMaMHPxJNia5OcndSe5K8sau/ZlJbkryte711PGVK0nqq8+I/zHgLVV1JnAucFmSM4GdwN6q2gzs7dYlSavEyMFfVYeq6vPd8n8C9wDrga3A7m633cBFPWuUJI3RWOb4k8wCzwduBU6rqkPdpoeA045xzI4k80nmFxYWxlGGJGkIvYM/yVOBvwJ+u6r+Y/G2qiqgjnZcVe2qqrmqmpuZmelbhiRpSL2CP8kTGYT+x6rqk13zN5Kc3m0/HTjcr0RJ0jj1+VRPgKuAe6rqPYs27QG2dcvbgBtGL0+SNG5rehz7QuBXgC8l+ULX9g7gCuC6JNuB+4GLe1UoSRqrkYO/qv4RyDE2nzfqeSVJk+Wdu5LUGINfkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNmVjwJ9mS5CtJ9ifZOanrSJKWZyLBn+Qk4APAy4EzgUuTnDmJa0mSlmdSI/5zgP1VdW9VfRf4OLB1QteSJC3Dmgmddz3w4KL1A8DPL94hyQ5gR7f6X0m+MqFaVtI64JvTLmKVsm+Ozb45thO+b3LlyIeuA35slAMnFfxLqqpdwK5pXX8SksxX1dy061iN7Jtjs2+Ozb45tq5vZkc5dlJTPQeBjYvWN3RtkqQpm1Tw/wuwOcmmJCcDlwB7JnQtSdIyTGSqp6oeS/KbwGeAk4Crq+quSVxrlTmhpq7GzL45Nvvm2OybYxu5b1JV4yxEkrTKeeeuJDXG4Jekxhj8y7TUoyiSvDnJ3UnuSLI3yUifsz1eDfuojiS/nKSSNPNRvWH6JsnF3e/PXUn+YqVrnJYh/l2dkeTmJLd3/7bOn0adKy3J1UkOJ7nzGNuT5P1dv92R5OyhTlxV/gz5w+CN6n8Ffhw4GfgicOYR+/wicEq3/HrgE9OuezX1T7ff04BbgH3A3LTrXi19A2wGbgdO7dafPe26V1Hf7AJe3y2fCdw37bpXqG9eBJwN3HmM7ecDnwYCnAvcOsx5HfEvz5KPoqiqm6vq0W51H4N7GFox7KM6/gC4EviflSxuyobpm18DPlBV3waoqsMrXOO0DNM3BTy9W34G8G8rWN/UVNUtwLceZ5etwEdrYB+wNsnpS53X4F+eoz2KYv3j7L+dwV/jVizZP91/RTdW1Y0rWdgqMMzvznOB5yb5pyT7kmxZseqma5i++T3g1UkOAJ8C3rAypa16y80kYIqPbDjRJXk1MAf8wrRrWS2SPAF4D/DaKZeyWq1hMN3zYgb/U7wlyc9U1cPTLGqVuBS4pqr+KMkLgD9PclZV/e+0CzseOeJfnqEeRZHkpcA7gQur6jsrVNtqsFT/PA04C/hskvsYzEnuaeQN3mF+dw4Ae6rqe1X1deCrDP4QnOiG6ZvtwHUAVfU54MkMHlLWupEej2PwL8+Sj6JI8nzgQwxCv5U52h943P6pqkeqal1Vzdbg4VL7GPTT/HTKXVHDPMbkbxiM9kmyjsHUz70rWOO0DNM3DwDnAST5KQbBv7CiVa5Oe4DXdJ/uORd4pKoOLXWQUz3LUMd4FEWS3wfmq2oP8G7gqcBfJgF4oKounFrRK2jI/mnSkH3zGeBlSe4Gvg/8TlX9+/SqXhlD9s1bgA8neRODN3pfW93HWk5kSa5lMBhY172/cTnwRICq+lMG73ecD+wHHgVeN9R5G+g7SdIiTvVIUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktSY/wOEJOh1ul8SBAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "draw_hist(yp[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAAPPklEQVR4nO3df6zddX3H8edLirpNJrhem64tu8zVbNXFQm4YxmVD2RRLYjHbCCRqNc1qDC66mWXV/aH7QVKzKYmJY9ZAqEbFbupoBptjjIW4DPSiWGkZ84pF2lV6FUQNkY363h/3S3ostz3n3nPPPeXD85GcnO/38/18z/d9P7n3db79nO/5NlWFJKktzxp3AZKkpWe4S1KDDHdJapDhLkkNMtwlqUErxl0AwMqVK2tycnLcZUjS08pdd931naqamG/bKRHuk5OTTE9Pj7sMSXpaSfLAibY5LSNJDTLcJalBhrskNchwl6QG9Q33JM9N8sUkX02yL8mfde3nJLkzyUySTyd5dtf+nG59pts+OeKfQZJ0nEHO3B8HXlVVLwM2AhcnuQB4P3B1Vf0S8Aiwteu/FXika7+66ydJWkZ9w73m/LBbPb17FPAq4O+79l3Apd3y5m6dbvtFSbJUBUuS+htozj3JaUnuBo4AtwDfAL5XVU90XQ4Ca7rlNcCDAN32R4Gfm+c1tyWZTjI9Ozs71A8hSfpJA4V7VR2tqo3AWuB84JeHPXBV7ayqqaqampiY9wtWkqRFWtA3VKvqe0luA14OnJlkRXd2vhY41HU7BKwDDiZZATwf+O4S1vwTJrffNKqX7uvAjkvGdmxJOplBrpaZSHJmt/xTwG8D9wK3Ab/bddsC3Ngt7+nW6bb/W/nfPUnSshrkzH01sCvJacy9Geyuqn9Msh+4IclfAl8Bru36Xwt8PMkM8DBw+QjqliSdRN9wr6q9wLnztN/P3Pz78e0/An5vSaqTJC2K31CVpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUoL7hnmRdktuS7E+yL8k7uvb3JTmU5O7usalnn3cnmUlyX5LXjPIHkCQ91YoB+jwBvKuqvpzkDOCuJLd0266uqr/u7ZxkA3A58BLg54F/TfLiqjq6lIVLkk6s75l7VR2uqi93yz8A7gXWnGSXzcANVfV4VX0TmAHOX4piJUmDWdCce5JJ4Fzgzq7p7Un2JrkuyVld2xrgwZ7dDjLPm0GSbUmmk0zPzs4uvHJJ0gkNHO5Jngd8BnhnVX0fuAZ4EbAROAx8YCEHrqqdVTVVVVMTExML2VWS1MdA4Z7kdOaC/RNV9VmAqnqoqo5W1Y+Bj3Js6uUQsK5n97VdmyRpmQxytUyAa4F7q+qDPe2re7q9HrinW94DXJ7kOUnOAdYDX1y6kiVJ/QxytcwrgDcCX0tyd9f2HuCKJBuBAg4AbwWoqn1JdgP7mbvS5kqvlJGk5dU33KvqC0Dm2XTzSfa5CrhqiLokSUPwG6qS1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIa1Dfck6xLcluS/Un2JXlH1/6CJLck+Xr3fFbXniQfSjKTZG+S80b9Q0iSftIgZ+5PAO+qqg3ABcCVSTYA24Fbq2o9cGu3DvBaYH332AZcs+RVS5JOqm+4V9Xhqvpyt/wD4F5gDbAZ2NV12wVc2i1vBj5Wc+4AzkyyeqkLlySd2ILm3JNMAucCdwKrqupwt+nbwKpueQ3wYM9uB7u2419rW5LpJNOzs7MLrVuSdBIDh3uS5wGfAd5ZVd/v3VZVBdRCDlxVO6tqqqqmJiYmFrKrJKmPgcI9yenMBfsnquqzXfNDT063dM9HuvZDwLqe3dd2bZKkZTLI1TIBrgXuraoP9mzaA2zplrcAN/a0v6m7auYC4NGe6RtJ0jJYMUCfVwBvBL6W5O6u7T3ADmB3kq3AA8Bl3babgU3ADPAY8JalLFiS1F/fcK+qLwA5weaL5ulfwJVD1iVJGoLfUJWkBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDVoxbgLkKRxm9x+09iOfWDHJSN53b5n7kmuS3IkyT09be9LcijJ3d1jU8+2dyeZSXJfkteMpGpJ0kkNMi1zPXDxPO1XV9XG7nEzQJINwOXAS7p9/ibJaUtVrCRpMH3DvapuBx4e8PU2AzdU1eNV9U1gBjh/iPokSYswzAeqb0+yt5u2OatrWwM82NPnYNcmSVpGiw33a4AXARuBw8AHFvoCSbYlmU4yPTs7u8gyJEnzWVS4V9VDVXW0qn4MfJRjUy+HgHU9Xdd2bfO9xs6qmqqqqYmJicWUIUk6gUWFe5LVPauvB568kmYPcHmS5yQ5B1gPfHG4EiVJC9X3OvcknwIuBFYmOQi8F7gwyUaggAPAWwGqal+S3cB+4Angyqo6OpLKJUkn1Dfcq+qKeZqvPUn/q4CrhilKkjQcbz8gSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIatGLcBTydTW6/aSzHPbDjkrEcV9LTh2fuktQgw12SGmS4S1KD+oZ7kuuSHElyT0/bC5LckuTr3fNZXXuSfCjJTJK9Sc4bZfGSpPkNcuZ+PXDxcW3bgVuraj1wa7cO8FpgfffYBlyzNGVKkhaib7hX1e3Aw8c1bwZ2dcu7gEt72j9Wc+4AzkyyeolqlSQNaLFz7quq6nC3/G1gVbe8Bniwp9/Bru0pkmxLMp1kenZ2dpFlSJLmM/QHqlVVQC1iv51VNVVVUxMTE8OWIUnqsdhwf+jJ6Zbu+UjXfghY19NvbdcmSVpGiw33PcCWbnkLcGNP+5u6q2YuAB7tmb6RJC2TvrcfSPIp4EJgZZKDwHuBHcDuJFuBB4DLuu43A5uAGeAx4C0jqFmS1EffcK+qK06w6aJ5+hZw5bBFSZKG4zdUJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBK8ZdgDSIye03je3YB3ZcMrZjS4s1VLgnOQD8ADgKPFFVU0leAHwamAQOAJdV1SPDlSlJWoilOHN/ZVV9p2d9O3BrVe1Isr1b/5MlOI46nsVK6mcUc+6bgV3d8i7g0hEcQ5J0EsOeuRfwL0kK+EhV7QRWVdXhbvu3gVXz7ZhkG7AN4Oyzzx6yDElLxX8ZtmHYcP/1qjqU5IXALUn+q3djVVUX/E/RvRHsBJiampq3j0494/zDlzS4oaZlqupQ93wE+BxwPvBQktUA3fORYYuUJC3MosM9yc8kOePJZeDVwD3AHmBL120LcOOwRUqSFmaYaZlVwOeSPPk6n6yqf07yJWB3kq3AA8Blw5cpSVqIRYd7Vd0PvGye9u8CFw1TlCRpOH5DVepjXB8ie+WIhuG9ZSSpQZ65S6coLzvVMDxzl6QGGe6S1CDDXZIaZLhLUoP8QFXSKcMPkZeOZ+6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkho0snBPcnGS+5LMJNk+quNIkp5qJOGe5DTgw8BrgQ3AFUk2jOJYkqSnGtWZ+/nATFXdX1X/C9wAbB7RsSRJx1kxotddAzzYs34Q+LXeDkm2Adu61R8muW9EtYzLSuA74y7iFOFYzHEcjnEsOnn/UGPxCyfaMKpw76uqdgI7x3X8UUsyXVVT467jVOBYzHEcjnEsjhnVWIxqWuYQsK5nfW3XJklaBqMK9y8B65Ock+TZwOXAnhEdS5J0nJFMy1TVE0neDnweOA24rqr2jeJYp7Bmp5wWwbGY4zgc41gcM5KxSFWN4nUlSWPkN1QlqUGGuyQ1yHAfUr/bLCT5oyT7k+xNcmuSE16X+nQ26O0mkvxOkkrS7GVwg4xFksu634t9ST653DUulwH+Ps5OcluSr3R/I5vGUeeoJbkuyZEk95xge5J8qBunvUnOG/qgVeVjkQ/mPiz+BvCLwLOBrwIbjuvzSuCnu+W3AZ8ed93jGIeu3xnA7cAdwNS46x7j78R64CvAWd36C8dd9xjHYifwtm55A3Bg3HWPaCx+AzgPuOcE2zcB/wQEuAC4c9hjeuY+nL63Waiq26rqsW71Duau+W/NoLeb+Avg/cCPlrO4ZTbIWPw+8OGqegSgqo4sc43LZZCxKOBnu+XnA/+zjPUtm6q6HXj4JF02Ax+rOXcAZyZZPcwxDffhzHebhTUn6b+VuXfn1vQdh+6fmeuq6qblLGwMBvmdeDHw4iT/keSOJBcvW3XLa5CxeB/whiQHgZuBP1ie0k45C82SvsZ2+4FnmiRvAKaA3xx3LcstybOADwJvHnMpp4oVzE3NXMjcv+RuT/KrVfW9cRY1JlcA11fVB5K8HPh4kpdW1Y/HXdjTnWfuwxnoNgtJfgv4U+B1VfX4MtW2nPqNwxnAS4F/T3KAuTnFPY1+qDrI78RBYE9V/V9VfRP4b+bCvjWDjMVWYDdAVf0n8Fzmbir2TLPkt2wx3IfT9zYLSc4FPsJcsLc6t3rScaiqR6tqZVVNVtUkc589vK6qpsdT7kgNcuuNf2DurJ0kK5mbprl/GWtcLoOMxbeAiwCS/Apz4T67rFWeGvYAb+qumrkAeLSqDg/zgk7LDKFOcJuFJH8OTFfVHuCvgOcBf5cE4FtV9bqxFT0CA47DM8KAY/F54NVJ9gNHgT+uqu+Or+rRGHAs3gV8NMkfMvfh6puru3ykJUk+xdwb+sru84X3AqcDVNXfMvd5wyZgBngMeMvQx2xwHCXpGc9pGUlqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGvT/Crz8cU2EFKgAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "draw_hist(yp[1])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD6CAYAAABamQdMAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAAQfElEQVR4nO3df4xlZX3H8fdHFtFWKtAdyXZ326F2jUUbFzKlGJsWoSouiYtpSyBRkJCuNdhoa4yr/UNtS4JplYTE0q6BshgVt/4oG8FaihhiU8BB15VdpK6wyG5XdpQfSoi0rN/+MYc67M7svTN37tz14f1Kbuac5zznnO88mf3MmeeeezZVhSSpLc8ZdQGSpMVnuEtSgwx3SWqQ4S5JDTLcJalBhrskNahnuCd5XpI7k3wzyY4kH+zar01yf5Jt3Wtt154kVybZlWR7klOH/D1Ikg6yrI8+TwJnVtXjSY4Gvprki922d1fVZw7q/3pgTff6HeCq7uucli9fXuPj4/MqXJKe7e66664fVNXYbNt6hntNf8rp8W716O51uE8+rQeu6/a7PclxSVZU1b65dhgfH2dycrJXKZKkGZI8MNe2vubckxyVZBuwH7i5qu7oNl3WTb1ckeSYrm0l8OCM3fd0bZKkJdJXuFfVgapaC6wCTkvycuC9wEuB3wZOAN4znxMn2ZBkMsnk1NTU/KqWJB3WvO6WqapHgVuBs6tqX017Evgn4LSu215g9YzdVnVtBx9rU1VNVNXE2NisU0aSpAXq526ZsSTHdcvPB14DfDvJiq4twLnA3d0uW4ELu7tmTgceO9x8uyRp8fVzt8wKYHOSo5j+ZbClqr6Q5MtJxoAA24A/7frfBKwDdgFPABcvetWSpMPq526Z7cAps7SfOUf/Ai4dvDRJ0kL5CVVJapDhLkkNMtwlqUH9vKEqSU0b33jjyM69+/JzhnJcr9wlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDWoZ7gneV6SO5N8M8mOJB/s2k9KckeSXUk+neS5Xfsx3fqubvv4kL8HSdJB+rlyfxI4s6peAawFzk5yOvAh4Iqq+g3gEeCSrv8lwCNd+xVdP0nSEuoZ7jXt8W716O5VwJnAZ7r2zcC53fL6bp1u+1lJslgFS5J662vOPclRSbYB+4Gbge8Cj1bVU12XPcDKbnkl8CBAt/0x4JdnOeaGJJNJJqempgb6JiRJz9RXuFfVgapaC6wCTgNeOuiJq2pTVU1U1cTY2Nigh5MkzTCvu2Wq6lHgVuCVwHFJlnWbVgF7u+W9wGqAbvsLgR8uRrGSpP70c7fMWJLjuuXnA68B7mE65P+o63YRcEO3vLVbp9v+5aqqRaxZktTDst5dWAFsTnIU078MtlTVF5LsBK5P8jfAN4Cru/5XAx9Psgt4GDh/CHVLkg6jZ7hX1XbglFna72N6/v3g9p8Af7wo1UmSFsRPqEpSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAb1DPckq5PcmmRnkh1J3tG1fyDJ3iTbute6Gfu8N8muJPcmed0wvwFJ0qGW9dHnKeBdVfX1JMcCdyW5udt2RVX93czOSU4GzgdeBvwK8O9JXlJVBxazcEnS3HpeuVfVvqr6erf8Y+AeYOVhdlkPXF9VT1bV/cAu4LTFKFaS1J95zbknGQdOAe7omt6eZHuSa5Ic37WtBB6csdseZvllkGRDkskkk1NTU/OvXJI0p77DPckLgM8C76yqHwFXAS8G1gL7gA/P58RVtamqJqpqYmxsbD67SpJ66CvckxzNdLB/oqo+B1BVD1XVgar6KfAxfjb1shdYPWP3VV2bJGmJ9HO3TICrgXuq6iMz2lfM6PZG4O5ueStwfpJjkpwErAHuXLySJUm99HO3zKuANwPfSrKta3sfcEGStUABu4G3AlTVjiRbgJ1M32lzqXfKSNLS6hnuVfVVILNsuukw+1wGXDZAXZKkAfgJVUlqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QG9Qz3JKuT3JpkZ5IdSd7RtZ+Q5OYk3+m+Ht+1J8mVSXYl2Z7k1GF/E5KkZ+rnyv0p4F1VdTJwOnBpkpOBjcAtVbUGuKVbB3g9sKZ7bQCuWvSqJUmHtaxXh6raB+zrln+c5B5gJbAeOKPrthn4CvCerv26qirg9iTHJVnRHUc/58Y33jiS8+6+/JyRnFf6eTWvOfck48ApwB3AiTMC+/vAid3ySuDBGbvt6dokSUuk73BP8gLgs8A7q+pHM7d1V+k1nxMn2ZBkMsnk1NTUfHaVJPXQV7gnOZrpYP9EVX2ua34oyYpu+wpgf9e+F1g9Y/dVXdszVNWmqpqoqomxsbGF1i9JmkU/d8sEuBq4p6o+MmPTVuCibvki4IYZ7Rd2d82cDjzmfLskLa2eb6gCrwLeDHwrybau7X3A5cCWJJcADwDnddtuAtYBu4AngIsXs2BJUm/93C3zVSBzbD5rlv4FXDpgXZKkAfgJVUlqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGLRt1AZq/8Y03jroESUe4nlfuSa5Jsj/J3TPaPpBkb5Jt3WvdjG3vTbIryb1JXjeswiVJc+tnWuZa4OxZ2q+oqrXd6yaAJCcD5wMv6/b5+yRHLVaxkqT+9Az3qroNeLjP460Hrq+qJ6vqfmAXcNoA9UmSFmCQN1TfnmR7N21zfNe2EnhwRp89XZskaQktNNyvAl4MrAX2AR+e7wGSbEgymWRyampqgWVIkmazoHCvqoeq6kBV/RT4GD+betkLrJ7RdVXXNtsxNlXVRFVNjI2NLaQMSdIcFhTuSVbMWH0j8PSdNFuB85Mck+QkYA1w52AlSpLmq+d97kk+BZwBLE+yB3g/cEaStUABu4G3AlTVjiRbgJ3AU8ClVXVgKJVLkubUM9yr6oJZmq8+TP/LgMsGKUqSNBgfPyBJDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhrUM9yTXJNkf5K7Z7SdkOTmJN/pvh7ftSfJlUl2Jdme5NRhFi9Jml0/V+7XAmcf1LYRuKWq1gC3dOsArwfWdK8NwFWLU6YkaT56hntV3QY8fFDzemBzt7wZOHdG+3U17XbguCQrFqlWSVKfFjrnfmJV7euWvw+c2C2vBB6c0W9P13aIJBuSTCaZnJqaWmAZkqTZDPyGalUVUAvYb1NVTVTVxNjY2KBlSJJmWGi4P/T0dEv3dX/XvhdYPaPfqq5NkrSEFhruW4GLuuWLgBtmtF/Y3TVzOvDYjOkbSdISWdarQ5JPAWcAy5PsAd4PXA5sSXIJ8ABwXtf9JmAdsAt4Arh4CDVLknroGe5VdcEcm86apW8Blw5alCRpMH5CVZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDWo562Q0pFgfOONIzv37svPGdm5pYXyyl2SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGDfQ89yS7gR8DB4CnqmoiyQnAp4FxYDdwXlU9MliZkqT5WIwr91dX1dqqmujWNwK3VNUa4JZuXZK0hIYxLbMe2NwtbwbOHcI5JEmHMWi4F/BvSe5KsqFrO7Gq9nXL3wdOnG3HJBuSTCaZnJqaGrAMSdJMg/4fqr9bVXuTvAi4Ocm3Z26sqkpSs+1YVZuATQATExOz9pEkLcxAV+5Vtbf7uh/4PHAa8FCSFQDd1/2DFilJmp8Fh3uSX0xy7NPLwGuBu4GtwEVdt4uAGwYtUpI0P4NMy5wIfD7J08f5ZFX9a5KvAVuSXAI8AJw3eJmSpPlYcLhX1X3AK2Zp/yFw1iBFSZIG4ydUJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYM+FfJZbXzjjaMuQZJm5ZW7JDXIK3eph1H9hbb78nNGcl61wXCXdMRwqnPxGO6SnsGAbYNz7pLUIMNdkhrktIx0hHJ6RIPwyl2SGmS4S1KDDHdJapDhLkkNGlq4Jzk7yb1JdiXZOKzzSJIONZS7ZZIcBXwUeA2wB/hakq1VtXOxz+UdBZJ0qGFduZ8G7Kqq+6rqf4DrgfVDOpck6SDDCveVwIMz1vd0bZKkJTCyDzEl2QBs6FYfT3LvqGqZp+XAD0ZdxBHGMZmd43Iox+Qg+dBAY/Jrc20YVrjvBVbPWF/Vtf2/qtoEbBrS+YcmyWRVTYy6jiOJYzI7x+VQjsmhhjUmw5qW+RqwJslJSZ4LnA9sHdK5JEkHGcqVe1U9leTtwJeAo4BrqmrHMM4lSTrU0Obcq+om4KZhHX+Efu6mkpaAYzI7x+VQjsmhhjImqaphHFeSNEI+fkCSGmS4z6HX4xOS/EWSnUm2J7klyZy3JLWi30dKJPnDJJWk+bsi+hmTJOd1Pys7knxyqWtcan382/nVJLcm+Ub372fdKOpcSkmuSbI/yd1zbE+SK7sx257k1IFPWlW+Dnox/Sbwd4FfB54LfBM4+aA+rwZ+oVt+G/DpUdc96jHp+h0L3AbcDkyMuu5RjwmwBvgGcHy3/qJR130EjMkm4G3d8snA7lHXvQTj8nvAqcDdc2xfB3wRCHA6cMeg5/TKfXY9H59QVbdW1RPd6u1M38vfsn4fKfHXwIeAnyxlcSPSz5j8CfDRqnoEoKr2L3GNS62fMSngl7rlFwL/vYT1jURV3QY8fJgu64HratrtwHFJVgxyTsN9dvN9fMIlTP/WbVnPMen+lFxdVc+Wp7n183PyEuAlSf4jye1Jzl6y6kajnzH5APCmJHuYvqPuz5amtCPaoj+yxf9DdUBJ3gRMAL8/6lpGKclzgI8AbxlxKUeaZUxPzZzB9F93tyX5rap6dJRFjdgFwLVV9eEkrwQ+nuTlVfXTURfWEq/cZ9fz8QkASf4A+EvgDVX15BLVNiq9xuRY4OXAV5LsZnrecGvjb6r283OyB9haVf9bVfcD/8V02LeqnzG5BNgCUFX/CTyP6WfOPJv1lTnzYbjPrufjE5KcAvwj08He+jwq9BiTqnqsqpZX1XhVjTP9PsQbqmpyNOUuiX4es/EvTF+1k2Q509M09y1hjUutnzH5HnAWQJLfZDrcp5a0yiPPVuDC7q6Z04HHqmrfIAd0WmYWNcfjE5L8FTBZVVuBvwVeAPxzEoDvVdUbRlb0kPU5Js8qfY7Jl4DXJtkJHADeXVU/HF3Vw9XnmLwL+FiSP2f6zdW3VHfLSKuSfIrpX/LLu/ca3g8cDVBV/8D0ew/rgF3AE8DFA5+z8TGVpGclp2UkqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDfo/Km92Kmk6CJQAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "draw_hist(yp[2])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAARGklEQVR4nO3dfYxld13H8ffHrgVBoIUdmrK7OCtsxYVgqEOtIRKwRkpBtgZCthFdtLoRKz5AAgWMJRpiASNCRJKV1i4GC01Fu1pQS2ltNLY4hdKyLQ9LKXTXlh3kSSEWCl//uAdzM53de+eeuXOnP96vZDPn/M6593x6Mv3s2d99OKkqJElt+b5ZB5AkrT3LXZIaZLlLUoMsd0lqkOUuSQ3aNOsAAJs3b675+flZx5CkB5Wbb775i1U1t9K2DVHu8/PzLC4uzjqGJD2oJPncsbY5LSNJDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ3aEJ9QlUaZv/DqmR37roufN7NjS5Pyyl2SGmS5S1KDLHdJapDlLkkNstwlqUEjyz3JpUmOJvn4svGXJ/lEkoNJ3jQ0/pokh5J8MslzphFaknR847wV8jLgz4B3fXcgybOBXcCPVdV9SR7bje8EdgNPBh4HfDDJaVX17bUOrtmY5VsSJY1v5JV7Vd0AfGnZ8MuAi6vqvm6fo934LuA9VXVfVX0WOAScsYZ5JUljmHTO/TTgp5LclORfkjy9G98C3D203+FuTJK0jib9hOom4NHAmcDTgSuS/PBqniDJXmAvwOMf//gJY0iSVjLplfth4H018GHgO8Bm4AiwbWi/rd3YA1TVvqpaqKqFubkVb94tSZrQpOX+d8CzAZKcBpwIfBE4AOxO8pAk24EdwIfXIKckaRVGTsskuRx4FrA5yWHgIuBS4NLu7ZHfBPZUVQEHk1wB3A7cD1zgO2Ukaf2NLPeqOu8Ym15yjP3fALyhTyhJUj9+QlWSGmS5S1KDLHdJapDlLkkNstwlqUGWuyQ1yHKXpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDbLcJalBlrskNchyl6QGWe6S1KCR5Z7k0iRHu7suLd/2yiSVZHO3niRvS3Ioya1JTp9GaEnS8Y1z5X4ZcPbywSTbgJ8FPj80/FwG903dAewF3tE/oiRptUaWe1XdAHxphU1vAV4F1NDYLuBdNXAjcFKSU9ckqSRpbBPNuSfZBRypqo8t27QFuHto/XA3ttJz7E2ymGRxaWlpkhiSpGNYdbkneRjwWuD3+xy4qvZV1UJVLczNzfV5KknSMpsmeMwTgO3Ax5IAbAU+kuQM4AiwbWjfrd2YJGkdrfrKvapuq6rHVtV8Vc0zmHo5varuBQ4Av9S9a+ZM4KtVdc/aRpYkjTLOWyEvB/4d+JEkh5Ocf5zd3w/cCRwC/gL4jTVJKUlalZHTMlV13ojt80PLBVzQP5YkqQ8/oSpJDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQZa7JDXIcpekBlnuktQgy12SGmS5S1KDLHdJatA4d2K6NMnRJB8fGntzkk8kuTXJ3yY5aWjba5IcSvLJJM+ZUm5J0nGMc+V+GXD2srFrgKdU1VOBTwGvAUiyE9gNPLl7zJ8nOWHN0kqSxjKy3KvqBuBLy8b+uaru71ZvBLZ2y7uA91TVfVX1WQb3Uj1jDfNKksawFnPuvwJ8oFveAtw9tO1wN/YASfYmWUyyuLS0tAYxJEnf1avck7wOuB9492ofW1X7qmqhqhbm5ub6xJAkLbNp0gcmeSnwfOCsqqpu+AiwbWi3rd2YJGkdTXTlnuRs4FXAC6rqG0ObDgC7kzwkyXZgB/Dh/jElSasx8so9yeXAs4DNSQ4DFzF4d8xDgGuSANxYVb9eVQeTXAHczmC65oKq+va0wkuSVjay3KvqvBWGLznO/m8A3tAnlI5v/sKrZx1B0gY38Zy79L1iVn+Z3nXx82ZyXLXBrx+QpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ0aWe5JLk1yNMnHh8YeneSaJJ/ufp7cjSfJ25IcSnJrktOnGV6StLJxrtwvA85eNnYhcG1V7QCu7dYBnsvg1no7gL3AO9YmpiRpNca5E9MNSeaXDe9icOs9gP3A9cCru/F3dTfMvjHJSUlOrap71iyx9D3Cm4Soj0nn3E8ZKux7gVO65S3A3UP7He7GJEnrqPcLqt1Veq32cUn2JllMsri0tNQ3hiRpyKTl/oUkpwJ0P49240eAbUP7be3GHqCq9lXVQlUtzM3NTRhDkrSSScv9ALCnW94DXDU0/kvdu2bOBL7qfLskrb+RL6gmuZzBi6ebkxwGLgIuBq5Icj7wOeDF3e7vB84BDgHfAH55CpklSSOM826Z846x6awV9i3ggr6hJEn9+AlVSWqQ5S5JDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQSO/W0bS95ZZ3QEKvAvUWvLKXZIaZLlLUoMsd0lqkOUuSQ3qVe5JfjfJwSQfT3J5kocm2Z7kpiSHkrw3yYlrFVaSNJ6Jyz3JFuC3gIWqegpwArAbeCPwlqp6IvBl4Py1CCpJGl/faZlNwA8k2QQ8DLgH+Gngym77fuDcnseQJK3SxOVeVUeAPwY+z6DUvwrcDHylqu7vdjsMbFnp8Un2JllMsri0tDRpDEnSCvpMy5wM7AK2A48DHg6cPe7jq2pfVS1U1cLc3NykMSRJK+gzLfMzwGeraqmqvgW8D3gGcFI3TQOwFTjSM6MkaZX6lPvngTOTPCxJgLOA24HrgBd1++wBruoXUZK0Wn3m3G9i8MLpR4DbuufaB7waeEWSQ8BjgEvWIKckaRV6fXFYVV0EXLRs+E7gjD7PK0nqx0+oSlKDLHdJapDlLkkNstwlqUGWuyQ1yHKXpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDbLcJalBlrskNchyl6QGWe6S1CDLXZIa1Kvck5yU5Mokn0hyR5KfTPLoJNck+XT38+S1CitJGk/fK/e3Av9YVU8Cfgy4A7gQuLaqdgDXduuSpHU0cbkneRTwTLp7pFbVN6vqK8AuYH+3237g3H4RJUmr1efKfTuwBPxlko8meWeShwOnVNU93T73Aqes9OAke5MsJllcWlrqEUOStFyfct8EnA68o6qeBnydZVMwVVVArfTgqtpXVQtVtTA3N9cjhiRpuT7lfhg4XFU3detXMij7LyQ5FaD7ebRfREnSak1c7lV1L3B3kh/phs4CbgcOAHu6sT3AVb0SSpJWbVPPx78ceHeSE4E7gV9m8BfGFUnOBz4HvLjnMSRJq9Sr3KvqFmBhhU1n9XneB4v5C6+edQRJWpGfUJWkBlnuktQgy12SGmS5S1KDLHdJapDlLkkNstwlqUGWuyQ1yHKXpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDbLcJalBvcs9yQndDbL/oVvfnuSmJIeSvLe7kYckaR2txZX7bwN3DK2/EXhLVT0R+DJw/hocQ5K0Cr3KPclW4HnAO7v1AD/N4GbZAPuBc/scQ5K0en2v3P8UeBXwnW79McBXqur+bv0wsGWlBybZm2QxyeLS0lLPGJKkYROXe5LnA0er6uZJHl9V+6pqoaoW5ubmJo0hSVpBnxtkPwN4QZJzgIcCjwTeCpyUZFN39b4VONI/piRpNSa+cq+q11TV1qqaB3YDH6qqXwCuA17U7bYHuKp3SknSqkzjfe6vBl6R5BCDOfhLpnAMSdJx9JmW+X9VdT1wfbd8J3DGWjyvJGkyfkJVkhpkuUtSgyx3SWqQ5S5JDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lq0Jp8t4wkrYX5C6+edYR1d9fFz5vK8z7oy/178ZdBkkZxWkaSGmS5S1KDLHdJalCfG2RvS3JdktuTHEzy2934o5Nck+TT3c+T1y6uJGkcfa7c7wdeWVU7gTOBC5LsBC4Erq2qHcC13bokaR31uUH2PVX1kW75v4E7gC3ALmB/t9t+4NyeGSVJq7Qmc+5J5oGnATcBp1TVPd2me4FTjvGYvUkWkywuLS2tRQxJUqd3uSf5QeBvgN+pqq8Nb6uqAmqlx1XVvqpaqKqFubm5vjEkSUN6lXuS72dQ7O+uqvd1w19Icmq3/VTgaL+IkqTV6vNumQCXAHdU1Z8MbToA7OmW9wBXTR5PkjSJPl8/8AzgF4HbktzSjb0WuBi4Isn5wOeAF/dKKElatYnLvar+FcgxNp816fNKkvrzE6qS1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQZa7JDXIcpekBlnuktQgy12SGmS5S1KDLHdJapDlLkkNstwlqUGWuyQ1yHKXpAZNrdyTnJ3kk0kOJblwWseRJD3QVMo9yQnA24HnAjuB85LsnMaxJEkPNK0r9zOAQ1V1Z1V9E3gPsGtKx5IkLdPnBtnHswW4e2j9MPATwzsk2Qvs7Vb/J8knV/H8m4Ev9kq4vsw7XeadLvNOUd7YK+8PHWvDtMp9pKraB+yb5LFJFqtqYY0jTY15p8u802Xe6ZpW3mlNyxwBtg2tb+3GJEnrYFrl/h/AjiTbk5wI7AYOTOlYkqRlpjItU1X3J/lN4J+AE4BLq+rgGh5ioumcGTLvdJl3usw7XVPJm6qaxvNKkmbIT6hKUoMsd0lq0IYu91FfYZDkmUk+kuT+JC+aRcZleUblfUWS25PcmuTaJMd8j+p6GCPvrye5LcktSf511p8yHvcrLZK8MEklmenb4cY4vy9NstSd31uS/Ooscg7lGXl+k7y4+x0+mOSv1zvjsiyjzu9bhs7tp5J8ZQYxh/OMyvv4JNcl+WjXEef0OmBVbcg/DF6I/Qzww8CJwMeAncv2mQeeCrwLeNGDIO+zgYd1yy8D3rvB8z5yaPkFwD9u5Lzdfo8AbgBuBBY2cl7gpcCfzSrjBHl3AB8FTu7WH7uR8y7b/+UM3tixYfMyeGH1Zd3yTuCuPsfcyFfuI7/CoKruqqpbge/MIuAy4+S9rqq+0a3eyOD9/7MyTt6vDa0+HJjlq+/jfqXFHwJvBP53PcOt4MH2FRzj5P014O1V9WWAqjq6zhmHrfb8ngdcvi7JVjZO3gIe2S0/CvjPPgfcyOW+0lcYbJlRlnGsNu/5wAemmuj4xsqb5IIknwHeBPzWOmVbyci8SU4HtlXV1esZ7BjG/X14YfdP8CuTbFth+3oZJ+9pwGlJ/i3JjUnOXrd0DzT2/2/d9Od24EPrkOtYxsn7euAlSQ4D72fwr42JbeRyb1aSlwALwJtnnWWUqnp7VT0BeDXwe7POcyxJvg/4E+CVs86yCn8PzFfVU4FrgP0zzjPKJgZTM89icCX8F0lOmmWgMe0Grqyqb886yAjnAZdV1VbgHOCvut/riWzkcn+wfYXBWHmT/AzwOuAFVXXfOmVbyWrP73uAc6cZaIRReR8BPAW4PsldwJnAgRm+qDry/FbVfw39DrwT+PF1yraScX4fDgMHqupbVfVZ4FMMyn4WVvP7u5vZTsnAeHnPB64AqKp/Bx7K4EvQJjOrFxjGeAFiE3Ang39OffcFiCcfY9/LmP0LqiPzAk9j8KLKjgfD+R3OCfwcsLiR8y7b/3pm+4LqOOf31KHlnwdu3OB5zwb2d8ubGUwzPGaj5u32exJwF90HNjf4+f0A8NJu+UcZzLlPnHtm/7FjnpBzGFwdfAZ4XTf2BwyuegGezuBq4uvAfwEHN3jeDwJfAG7p/hzY4HnfChzssl53vDLdCHmX7TvTch/z/P5Rd34/1p3fJ23wvGEw9XU7cBuweyPn7dZfD1w8y5yrOL87gX/rfh9uAX62z/H8+gFJatBGnnOXJE3IcpekBlnuktQgy12SGmS5S1KDLHdJapDlLkkN+j9slSDV7zxgkgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "draw_hist(yp[3])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAANsElEQVR4nO3df6zd9V3H8ecL6macTGDtmqatu3PpMjujgDdIosmYmAklW2EzpCRzZUGrC/NHnInVmYzMEDuNW1xEkm4jdIuD4XShBpxihZAtVneR8nMyCiuhXaF3v2eIU/DtH/dbdyi33HPvufecL589H8nJ/X4/3++559XvPX3lez/nnO9NVSFJasspkw4gSVp+lrskNchyl6QGWe6S1CDLXZIatGrSAQBWr15dU1NTk44hSS8qd99991eras1823pR7lNTU8zMzEw6hiS9qCR5/GTbnJaRpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDbLcJalBlrskNchyl6QG9eITqi9WUztvnXSEsTu06+JJR5A0BM/cJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQZa7JDXIcpekBlnuktQgy12SGmS5S1KDLHdJapDlLkkNstwlqUGWuyQ1yHKXpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDbLcJalBlrskNWjBck+yMckdSR5K8mCS3+rGz0xye5JHuq9ndONJ8uEkB5Pcl+Sclf5HSJKea5gz92eA91TVZuA84Kokm4GdwL6q2gTs69YBLgI2dbcdwHXLnlqS9IIWLPeqOlpV/94tfwf4IrAe2Ars6XbbA1zSLW8FPl5z9gOnJ1m33MElSSe3qDn3JFPA2cC/Amur6mi36Ulgbbe8Hnhi4G6HuzFJ0pgMXe5Jfhj4G+C3q+rbg9uqqoBazAMn2ZFkJsnM7OzsYu4qSVrAUOWe5AeYK/a/qqq/7YafOj7d0n091o0fATYO3H1DN/YcVbW7qqaranrNmjVLzS9Jmscw75YJ8DHgi1X1wYFNe4Ht3fJ24JaB8Xd075o5D/jWwPSNJGkMVg2xz88Cvwzcn+RAN/YHwC7g5iRXAo8Dl3XbbgO2AAeBp4F3LmdgSdLCFiz3qvockJNsvmCe/Qu4asRckqQR+AlVSWqQ5S5JDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lq0DCXH5D+39TOWyfyuId2XTyRx5VerDxzl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQZa7JDXIcpekBlnuktQgy12SGmS5S1KDLHdJapDlLkkNstwlqUGWuyQ1yHKXpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDbLcJalBlrskNchyl6QGLVjuSa5PcizJAwNjVyc5kuRAd9sysO33kxxM8nCSX1yp4JKkkxvmzP0G4MJ5xj9UVWd1t9sAkmwGtgGv7+7zl0lOXa6wkqThLFjuVXUX8PUhv99W4Kaq+m5VfRk4CJw7Qj5J0hKMMuf+7iT3ddM2Z3Rj64EnBvY53I1JksZoqeV+HfAa4CzgKPBni/0GSXYkmUkyMzs7u8QYkqT5LKncq+qpqnq2qv4X+Ajfm3o5Amwc2HVDNzbf99hdVdNVNb1mzZqlxJAkncSSyj3JuoHVS4Hj76TZC2xL8tIkrwY2Af82WkRJ0mKtWmiHJDcC5wOrkxwG3gecn+QsoIBDwK8BVNWDSW4GHgKeAa6qqmdXJLkk6aQWLPequnye4Y+9wP7XANeMEkqSNBo/oSpJDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQZa7JDXIcpekBlnuktQgy12SGmS5S1KDLHdJapDlLkkNstwlqUGWuyQ1yHKXpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQZa7JDXIcpekBi1Y7kmuT3IsyQMDY2cmuT3JI93XM7rxJPlwkoNJ7ktyzkqGlyTNb5gz9xuAC08Y2wnsq6pNwL5uHeAiYFN32wFctzwxJUmLsWC5V9VdwNdPGN4K7OmW9wCXDIx/vObsB05Psm6ZskqShrTUOfe1VXW0W34SWNstrweeGNjvcDf2PEl2JJlJMjM7O7vEGJKk+Yz8gmpVFVBLuN/uqpququk1a9aMGkOSNGCp5f7U8emW7uuxbvwIsHFgvw3dmCRpjJZa7nuB7d3yduCWgfF3dO+aOQ/41sD0jSRpTFYttEOSG4HzgdVJDgPvA3YBNye5EngcuKzb/TZgC3AQeBp45wpkliQtYMFyr6rLT7Lpgnn2LeCqUUNJkkbjJ1QlqUGWuyQ1yHKXpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDbLcJalBlrskNWjByw/03dTOWycdQZJ6xzN3SWqQ5S5JDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQZa7JDXIcpekBlnuktQgy12SGmS5S1KDLHdJapDlLkkNstwlqUGWuyQ1yHKXpAZZ7pLUIMtdkhpkuUtSg1aNcuckh4DvAM8Cz1TVdJIzgU8BU8Ah4LKq+sZoMSVJi7EcZ+5vrKqzqmq6W98J7KuqTcC+bl2SNEYrMS2zFdjTLe8BLlmBx5AkvYBRy72Af0xyd5Id3djaqjraLT8JrJ3vjkl2JJlJMjM7OztiDEnSoJHm3IGfq6ojSV4J3J7kPwY3VlUlqfnuWFW7gd0A09PT8+4jHTe189aJPfahXRdP7LGlpRqp3KvqSPf1WJLPAOcCTyVZV1VHk6wDji1DTklaMS2ePCx5WibJy5KcdnwZeBPwALAX2N7tth24ZdSQkqTFGeXMfS3wmSTHv88nq+qzSb4A3JzkSuBx4LLRY0qSFmPJ5V5VjwE/Nc/414ALRgklSRqNn1CVpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQZa7JDXIcpekBo30B7Kl7weT+uPJK/WHk/X9wTN3SWqQ5S5JDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQV5bRtJzTOpaOlpelrvUU5asRuG0jCQ1yHKXpAZZ7pLUIMtdkhq0YuWe5MIkDyc5mGTnSj2OJOn5VqTck5wKXAtcBGwGLk+yeSUeS5L0fCt15n4ucLCqHquq/wZuArau0GNJkk6wUu9zXw88MbB+GPiZwR2S7AB2dKv/meThRT7GauCrS044Hn3P2Pd8YMbl0Pd80P+MK5YvHxjp7q862YaJfYipqnYDu5d6/yQzVTW9jJGWXd8z9j0fmHE59D0f9D9j3/PNZ6WmZY4AGwfWN3RjkqQxWKly/wKwKcmrk7wE2AbsXaHHkiSdYEWmZarqmSTvBv4BOBW4vqoeXOaHWfKUzhj1PWPf84EZl0Pf80H/M/Y93/OkqiadQZK0zPyEqiQ1yHKXpAb1vtwXuoxBkt9J8lCS+5LsS3LS931OMOOvJ7k/yYEknxv3p3WHvRREkrclqSRjf8vXEMfwiiSz3TE8kORX+pSv2+ey7rn4YJJPjjPfMBmTfGjg+H0pyTd7lu9Hk9yR5J7u//OWceYbMuOrup65L8mdSTaMO+PQqqq3N+ZejH0U+DHgJcC9wOYT9nkj8EPd8ruAT/Uw48sHlt8CfLZP+br9TgPuAvYD0z08hlcAf9Hj5+Em4B7gjG79lX3LeML+v8HcGx16k4+5Fy3f1S1vBg717RgCfw1s75Z/HvjEJJ6Tw9z6fua+4GUMquqOqnq6W93P3Hvq+5bx2wOrLwPG+Sr2sJeC+CPgA8B/jTHbcX2/XMUw+X4VuLaqvgFQVcd6mHHQ5cCNY0k2Z5h8Bby8W/4R4CtjzAfDZdwM/HO3fMc823uj7+U+32UM1r/A/lcCf7+iiZ5vqIxJrkryKPAnwG+OKRsMkS/JOcDGqprU33Ub9uf8tu7X4U8n2TjP9pUyTL7XAq9N8vkk+5NcOLZ0c4b+v9JNXb6a75XUOAyT72rg7UkOA7cx99vFOA2T8V7grd3ypcBpSV4xhmyL1vdyH1qStwPTwJ9OOst8quraqnoN8HvAH046z3FJTgE+CLxn0lkW8HfAVFX9JHA7sGfCeU60irmpmfOZOyv+SJLTJxnoBWwDPl1Vz046yAkuB26oqg3AFuAT3fOzT34XeEOSe4A3MPfJ+74dR6D/5T7UZQyS/ALwXuAtVfXdMWU7brGXWrgJuGQlA51goXynAT8B3JnkEHAesHfML6oueAyr6msDP9uPAj89pmww3M/4MLC3qv6nqr4MfIm5sh+XxTwPtzHeKRkYLt+VwM0AVfUvwA8yd8GucRnmefiVqnprVZ3NXOdQVd8cW8LFmPSk/wIvcKwCHmPuV8jjL3C8/oR9zmbuRZBNPc64aWD5zcBMn/KdsP+djP8F1WGO4bqB5UuB/T3LdyGwp1tezdyv96/oU8Zuv9cBh+g+wNinfMxNqV7RLf84c3PuY8s5ZMbVwCnd8jXA+8d5HBf175l0gCEO+BbmzoIeBd7bjb2fubN0gH8CngIOdLe9Pcz458CDXb47XqhcJ5HvhH3HXu5DHsM/7o7hvd0xfF3P8oW56a2HgPuBbX07ht361cCucWcb8hhuBj7f/YwPAG/qYcZfAh7p9vko8NJJHMthbl5+QJIa1Pc5d0nSEljuktQgy12SGmS5S1KDLHdJapDlLkkNstwlqUH/B8/g61Hxq30HAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "draw_hist(yp[4])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAD4CAYAAAANbUbJAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAAQZklEQVR4nO3df4xldX3G8fcjVJtaKOiOhAA6QBYtGrvYKTVpVfxRi9CCqKFsqgVLXTHSptGmRW0qsTFdf6Cx0WIWJYARBKVEGtBKKUo00jrIuoKKsjjExXV3BH9VLRX49I85G6/j7M7duffOvcv3/Uomc873fu85T+7uPjl77rnnpqqQJD2yPWrcASRJo2fZS1IDLHtJaoBlL0kNsOwlqQH7jzsAwJo1a2p6enrcMSRpn3Lrrbd+t6qm+pk7EWU/PT3N7OzsuGNI0j4lyT39zvU0jiQ1wLKXpAZY9pLUAMtekhpg2UtSAyx7SWqAZS9JDbDsJakBlr0kNWAiPkE7qOnzrhvLfuc2njyW/UrS3vLIXpIaYNlLUgMse0lqgGUvSQ2w7CWpAZa9JDXAspekBlj2ktQAy16SGmDZS1IDli37JBcn2Znk9p6xK5Ns7n7mkmzuxqeT/LTnsfePMLskqU/93BvnEuC9wGW7BqrqT3YtJ7kA+EHP/K1VtW5I+SRJQ7Bs2VfVzUmml3osSYDTgecNOZckaYgGPWf/LGBHVX2jZ+zIJLcl+UySZ+3uiUk2JJlNMjs/Pz9gDEnSngxa9uuBK3rWtwNPrKrjgNcBlyc5cKknVtWmqpqpqpmpqakBY0iS9mTFZZ9kf+AlwJW7xqrqgaq6r1u+FdgKHDNoSEnSYAY5sn8B8LWq2rZrIMlUkv265aOAtcDdg0WUJA2qn0svrwA+Dzw5ybYkZ3cPncEvnsIBeDawpbsU82PAOVV1/xDzSpJWoJ+rcdbvZvysJcauBq4ePJYkaZj8BK0kNcCyl6QGWPaS1IB+bpeg3Zg+77qx7Hdu48lj2a+kfZdH9pLUAMtekhpg2UtSAyx7SWqAZS9JDbDsJakBlr0kNcCyl6QGWPaS1ADLXpIaYNlLUgMse0lqgGUvSQ2w7CWpAf18B+3FSXYmub1n7Pwk9ybZ3P2c1PPYG5LcleTOJH84quCSpP71c2R/CXDiEuPvrqp13c/1AEmOZeGLyJ/aPedfkuw3rLCSpJVZtuyr6mbg/j63dyrwkap6oKq+CdwFHD9APknSEAxyzv7cJFu60zwHd2OHAd/qmbOtG5MkjdFKy/5C4GhgHbAduGBvN5BkQ5LZJLPz8/MrjCFJ6seKyr6qdlTVQ1X1MHARPz9Vcy9wRM/Uw7uxpbaxqapmqmpmampqJTEkSX1aUdknObRn9TRg15U61wJnJHlMkiOBtcB/DxZRkjSo/ZebkOQK4ARgTZJtwJuBE5KsAwqYA14NUFV3JLkK+ArwIPDaqnpoJMklSX1btuyrav0Swx/cw/y3Am8dJJQkabj8BK0kNcCyl6QGWPaS1ADLXpIaYNlLUgMse0lqgGUvSQ2w7CWpAZa9JDXAspekBlj2ktQAy16SGmDZS1IDLHtJaoBlL0kNsOwlqQGWvSQ1wLKXpAZY9pLUgGXLPsnFSXYmub1n7B1JvpZkS5JrkhzUjU8n+WmSzd3P+0eYXZLUp36O7C8BTlw0dgPwtKp6OvB14A09j22tqnXdzznDiSlJGsSyZV9VNwP3Lxr7VFU92K3eAhw+gmySpCEZxjn7Pwc+0bN+ZJLbknwmybN296QkG5LMJpmdn58fQgxJ0u4MVPZJ3gQ8CHy4G9oOPLGqjgNeB1ye5MClnltVm6pqpqpmpqamBokhSVrGiss+yVnAHwF/WlUFUFUPVNV93fKtwFbgmCHklCQNYEVln+RE4G+BU6rqJz3jU0n265aPAtYCdw8jqCRp5fZfbkKSK4ATgDVJtgFvZuHqm8cANyQBuKW78ubZwFuS/Ax4GDinqu5fcsOSpFWzbNlX1folhj+4m7lXA1cPGkqTa/q868ay37mNJ49lv9IjhZ+glaQGWPaS1ADLXpIaYNlLUgMse0lqgGUvSQ2w7CWpAZa9JDXAspekBlj2ktQAy16SGmDZS1IDLHtJaoBlL0kNsOwlqQHL3s9ek2dc95SXtO/yyF6SGmDZS1ID+ir7JBcn2Znk9p6xxyW5Ick3ut8Hd+NJ8s9J7kqyJckzRhVektSffo/sLwFOXDR2HnBjVa0FbuzWAV4ErO1+NgAXDh5TkjSIvsq+qm4G7l80fCpwabd8KfDinvHLasEtwEFJDh1CVknSCg1yzv6QqtreLX8HOKRbPgz4Vs+8bd2YJGlMhvIGbVUVUHvznCQbkswmmZ2fnx9GDEnSbgxS9jt2nZ7pfu/sxu8FjuiZd3g39guqalNVzVTVzNTU1AAxJEnLGaTsrwXO7JbPBD7eM/5n3VU5zwR+0HO6R5I0Bn19gjbJFcAJwJok24A3AxuBq5KcDdwDnN5Nvx44CbgL+AnwyiFnliTtpb7KvqrW7+ah5y8xt4DXDhJKWmyct4iY23jy2PYtDYufoJWkBlj2ktQAy16SGmDZS1IDLHtJaoBlL0kNsOwlqQGWvSQ1wLKXpAZY9pLUAMtekhpg2UtSAyx7SWqAZS9JDbDsJakBlr0kNcCyl6QGWPaS1ADLXpIa0Nd30C4lyZOBK3uGjgL+ATgIeBUw342/saquX+l+JEmDW3HZV9WdwDqAJPsB9wLXAK8E3l1V7xxGQEnS4IZ1Guf5wNaqumdI25MkDdGwyv4M4Iqe9XOTbElycZKDl3pCkg1JZpPMzs/PLzVFkjQkA5d9kkcDpwAf7YYuBI5m4RTPduCCpZ5XVZuqaqaqZqampgaNIUnag2Ec2b8I+GJV7QCoqh1V9VBVPQxcBBw/hH1IkgYwjLJfT88pnCSH9jx2GnD7EPYhSRrAiq/GAUjyWOAPgFf3DL89yTqggLlFj0mSxmCgsq+qHwOPXzT2ioESSZKGzk/QSlIDLHtJaoBlL0kNGOicvaTRmT7vurHte27jyWPbt0bDI3tJaoBlL0kNsOwlqQGWvSQ1wLKXpAZY9pLUAMtekhpg2UtSAyx7SWqAZS9JDfB2CdIyxnnbAmlYPLKXpAZY9pLUAMtekhpg2UtSAwZ+gzbJHPAj4CHgwaqaSfI44EpgmoUvHT+9qr436L4kSSszrCP751bVuqqa6dbPA26sqrXAjd26JGlMRnUa51Tg0m75UuDFI9qPJKkPwyj7Aj6V5NYkG7qxQ6pqe7f8HeCQxU9KsiHJbJLZ+fn5IcSQJO3OMD5U9ftVdW+SJwA3JPla74NVVUlq8ZOqahOwCWBmZuaXHpckDc/AR/ZVdW/3eydwDXA8sCPJoQDd752D7keStHIDlX2SxyY5YNcy8ELgduBa4Mxu2pnAxwfZjyRpMIOexjkEuCbJrm1dXlWfTPIF4KokZwP3AKcPuB9J0gAGKvuquhv4rSXG7wOeP8i2JUnD4ydoJakBlr0kNcCyl6QGWPaS1ADLXpIaYNlLUgP8DlpJE2Nc3/c7t/Hksex3NXlkL0kNsOwlqQGWvSQ1wLKXpAZY9pLUAMtekhpg2UtSAyx7SWqAZS9JDbDsJakBlr0kNWDF98ZJcgRwGQvfQ1vApqp6T5LzgVcB893UN1bV9YMGlbR6xnWPGo3OIDdCexB4fVV9MckBwK1Jbugee3dVvXPweJKkYVhx2VfVdmB7t/yjJF8FDhtWMEnS8AzlnH2SaeA44L+6oXOTbElycZKDd/OcDUlmk8zOz88vNUWSNCQDl32SXweuBv66qn4IXAgcDaxj4cj/gqWeV1WbqmqmqmampqYGjSFJ2oOByj7Jr7BQ9B+uqn8FqKodVfVQVT0MXAQcP3hMSdIgVlz2SQJ8EPhqVb2rZ/zQnmmnAbevPJ4kaRgGuRrn94BXAF9OsrkbeyOwPsk6Fi7HnANePcA+JElDMMjVOJ8FssRDXlMvSRPGT9BKUgMse0lqgGUvSQ2w7CWpAYNcjSNJjwjjvPHb3MaTV2U/HtlLUgMse0lqgGUvSQ2w7CWpAZa9JDXAspekBlj2ktQAy16SGmDZS1IDLHtJaoBlL0kNsOwlqQGWvSQ1wLKXpAaMrOyTnJjkziR3JTlvVPuRJC1vJGWfZD/gfcCLgGOB9UmOHcW+JEnLG9WR/fHAXVV1d1X9H/AR4NQR7UuStIxRfVPVYcC3eta3Ab/bOyHJBmBDt/o/Se4c4v7XAN8d4vaGaVKzTWouMNtKmW1lVjVb3rZX0xdne1K/Txzb1xJW1SZg0yi2nWS2qmZGse1BTWq2Sc0FZlsps63MIzXbqE7j3Asc0bN+eDcmSRqDUZX9F4C1SY5M8mjgDODaEe1LkrSMkZzGqaoHk5wL/DuwH3BxVd0xin3txkhODw3JpGab1FxgtpUy28o8IrOlqoYZRJI0gfwErSQ1wLKXpAbss2W/3O0Ykjw7yReTPJjkZROW7XVJvpJkS5Ibk/R9rewqZDsnyZeTbE7y2dX85HO/t9hI8tIklWTVLo/r43U7K8l897ptTvIXk5Ktm3N693fujiSXT0q2JO/uec2+nuT7E5TtiUluSnJb92/1pAnK9qSuO7Yk+XSSw5fdaFXtcz8svOm7FTgKeDTwJeDYRXOmgacDlwEvm7BszwV+rVt+DXDlBGU7sGf5FOCTk5Ktm3cAcDNwCzAzKdmAs4D3rtbfs73Mtha4DTi4W3/CpGRbNP8vWbiYYyKysfBm6Gu65WOBuQnK9lHgzG75ecCHltvuvnpkv+ztGKpqrqq2AA9PYLabquon3eotLHwOYVKy/bBn9bHAar2D3+8tNv4ReBvwv6uUa2+yjUM/2V4FvK+qvgdQVTsnKFuv9cAVq5Ksv2wFHNgt/wbw7QnKdizwn93yTUs8/kv21bJf6nYMh40py2J7m+1s4BMjTfRzfWVL8tokW4G3A381KdmSPAM4oqquW6VMu/T7Z/rS7r/VH0tyxBKPj0I/2Y4BjknyuSS3JDlxgrIBC6clgCP5eYGNWj/ZzgdenmQbcD0L//NYDf1k+xLwkm75NOCAJI/f00b31bJ/REjycmAGeMe4s/SqqvdV1dHA3wF/P+48AEkeBbwLeP24s+zGvwHTVfV04Abg0jHn6bU/C6dyTmDh6PmiJAeNM9ASzgA+VlUPjTtIj/XAJVV1OHAS8KHu7+Ek+BvgOUluA57Dwh0K9vjaTUrwvTXJt2PoK1uSFwBvAk6pqgcmKVuPjwAvHmWgHstlOwB4GvDpJHPAM4FrV+lN2mVft6q6r+fP8QPAb69Crr6ysXBkeG1V/ayqvgl8nYXyn4Rsu5zB6p3Cgf6ynQ1cBVBVnwd+lYUbkY09W1V9u6peUlXHsdAjVNX397jV1XjDYQRvYOwP3M3Cf/t2vYHx1N3MvYTVfYN22WzAcSy8AbN20l633kzAHwOzk5Jt0fxPs3pv0Pbzuh3as3wacMsEZTsRuLRbXsPCKYLHT0K2bt5TgDm6D3lO0Ov2CeCsbvk3WThnP/KMfWZbAzyqW34r8JZlt7taL+4IXpCTWDhC2Qq8qRt7CwtHygC/w8IRzY+B+4A7JijbfwA7gM3dz7UTlO09wB1drpv2VLirnW3R3FUr+z5ft3/qXrcvda/bUyYoW1g4BfYV4MvAGZOSrVs/H9i4Wpn24nU7Fvhc92e6GXjhBGV7GfCNbs4HgMcst01vlyBJDdhXz9lLkvaCZS9JDbDsJakBlr0kNcCyl6QGWPaS1ADLXpIa8P9Ft+nig/k+FQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "draw_hist(yp[5])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAANiklEQVR4nO3dfaze5V3H8fdndIDGDRg9EtJWD8tYtJlxkIaxLNEJangwlERGIM7h0thsMjODiVb3h49/wB8OJSHTRsjKogNEI82GMchDiItlHoTxmG1nCNLK6BkP1YUwh/v6x32xHGpP77s990N79f1KTs71u37Xff++V+/TT3/nun/3r6kqJEl9ecusC5AkjZ/hLkkdMtwlqUOGuyR1yHCXpA6tmXUBAGvXrq35+flZlyFJR5WHHnroW1U1d6B9R0S4z8/Ps7CwMOsyJOmokuTZlfa5LCNJHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR06Ij6huhrz2744s2M/c+3FMzu2JB2MZ+6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOjRzuSY5L8nCSL7TtM5I8mGQxyW1Jjm/9J7TtxbZ/fkK1S5JWcChn7p8Enlq2fR1wfVW9C3gZ2NL6twAvt/7r2zhJ0hSNFO5J1gMXA3/ZtgOcB9zRhuwALm3tzW2btv/8Nl6SNCWjnrn/KfBbwPfa9qnAK1X1etveDaxr7XXAcwBt/742/k2SbE2ykGRhaWnp8KqXJB3Q0HBP8gvA3qp6aJwHrqrtVbWpqjbNzc2N86kl6Zi3ZoQxHwAuSXIRcCLwduDPgJOTrGln5+uBPW38HmADsDvJGuAk4MWxVy5JWtHQM/eq+p2qWl9V88AVwL1V9UvAfcBlbdhVwJ2tvbNt0/bfW1U11qolSQe1muvcfxu4JskigzX1m1r/TcCprf8aYNvqSpQkHapRlmW+r6ruB+5v7aeBcw4w5jXgQ2OoTZJ0mPyEqiR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1KGh4Z7kxCRfTvKVJE8k+YPWf0aSB5MsJrktyfGt/4S2vdj2z094DpKk/Yxy5v4d4Lyq+kngvcAFSc4FrgOur6p3AS8DW9r4LcDLrf/6Nk6SNEVDw70Gvt0239q+CjgPuKP17wAube3NbZu2//wkGVfBkqThRlpzT3JckkeAvcDdwDeAV6rq9TZkN7CutdcBzwG0/fuAU8dYsyRpiJHCvar+t6reC6wHzgF+bLUHTrI1yUKShaWlpdU+nSRpmUO6WqaqXgHuA94PnJxkTdu1HtjT2nuADQBt/0nAiwd4ru1VtamqNs3NzR1e9ZKkAxrlapm5JCe39g8APwc8xSDkL2vDrgLubO2dbZu2/96qqjHWLEkaYs3wIZwO7EhyHIN/DG6vqi8keRK4NckfAw8DN7XxNwGfS7IIvARcMYG6JUkHMTTcq+pR4KwD9D/NYP19//7XgA+NpTpJ0mHxE6qS1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6tDQcE+yIcl9SZ5M8kSST7b+dyS5O8nX2/dTWn+S3JBkMcmjSc6e9CQkSW82ypn768BvVtVG4Fzg6iQbgW3APVV1JnBP2wa4EDizfW0FPjP2qiVJBzU03Kvq+ar6t9b+b+ApYB2wGdjRhu0ALm3tzcAtNbALODnJ6eMuXJK0skNac08yD5wFPAicVlXPt13fBE5r7XXAc8setrv17f9cW5MsJFlYWlo61LolSQcxcrgn+SHgb4HfqKr/Wr6vqgqoQzlwVW2vqk1VtWlubu5QHipJGmKkcE/yVgbB/ldV9Xet+4U3llva972tfw+wYdnD17c+SdKUjHK1TICbgKeq6tPLdu0Ermrtq4A7l/V/pF01cy6wb9nyjSRpCtaMMOYDwC8DjyV5pPX9LnAtcHuSLcCzwOVt313ARcAi8Crw0XEWLEkabmi4V9U/A1lh9/kHGF/A1ausS5K0Cn5CVZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdGhruSW5OsjfJ48v63pHk7iRfb99Paf1JckOSxSSPJjl7ksVLkg5slDP3zwIX7Ne3Dbinqs4E7mnbABcCZ7avrcBnxlOmJOlQDA33qnoAeGm/7s3AjtbeAVy6rP+WGtgFnJzk9DHVKkka0eGuuZ9WVc+39jeB01p7HfDcsnG7W9//k2RrkoUkC0tLS4dZhiTpQFb9hmpVFVCH8bjtVbWpqjbNzc2ttgxJ0jKHG+4vvLHc0r7vbf17gA3Lxq1vfZKkKVpzmI/bCVwFXNu+37ms/xNJbgXeB+xbtnwjSUek+W1fnNmxn7n24ok879BwT/J54IPA2iS7gd9jEOq3J9kCPAtc3obfBVwELAKvAh+dQM2SpCGGhntVXbnCrvMPMLaAq1dblCRpdfyEqiR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDh3v7Ac1Qjx+VljRenrlLUoc8c1+FWZ5BS9LBeOYuSR0y3CWpQ4a7JHXIcJekDhnuktQhr5aR9CZ+jqIPnrlLUocMd0nqkOEuSR1yzV2HZFbrsa7FSofGM3dJ6pDhLkkdMtwlqUOuueuocCyu9XvXUa2GZ+6S1CHDXZI6ZLhLUocMd0nqkG+oSgfhm5o6WnnmLkkdMtwlqUOGuyR1yDV3SUcM3+MYH8/cJalDhrskdWgi4Z7kgiRfTbKYZNskjiFJWtnYwz3JccCNwIXARuDKJBvHfRxJ0somceZ+DrBYVU9X1f8AtwKbJ3AcSdIKJnG1zDrguWXbu4H37T8oyVZga9v8dpKvHuQ51wLfGluFR5djde7O+9hyrM6bXLequf/oSjtmdilkVW0Hto8yNslCVW2acElHpGN17s772HKszhsmN/dJLMvsATYs217f+iRJUzKJcP9X4MwkZyQ5HrgC2DmB40iSVjD2ZZmqej3JJ4B/BI4Dbq6qJ1b5tCMt33TqWJ278z62HKvzhgnNPVU1ieeVJM2Qn1CVpA4Z7pLUoSMq3IfdtiDJCUlua/sfTDI/gzLHboR5X5PkySSPJrknyYrXth5tRr1VRZJfTFJJurhcbpR5J7m8ve5PJPnradc4CSP8rP9IkvuSPNx+3i+aRZ3jluTmJHuTPL7C/iS5of25PJrk7FUftKqOiC8Gb75+A3gncDzwFWDjfmN+Dfjz1r4CuG3WdU9p3j8D/GBrf7yHeY869zbubcADwC5g06zrntJrfibwMHBK2/7hWdc9pXlvBz7e2huBZ2Zd95jm/lPA2cDjK+y/CPgHIMC5wIOrPeaRdOY+ym0LNgM7WvsO4PwkmWKNkzB03lV1X1W92jZ3MfjsQA9GvVXFHwHXAa9Ns7gJGmXevwrcWFUvA1TV3inXOAmjzLuAt7f2ScB/TrG+iamqB4CXDjJkM3BLDewCTk5y+mqOeSSF+4FuW7BupTFV9TqwDzh1KtVNzijzXm4Lg3/hezB07u3X0w1V1dP/4jDKa/5u4N1JvpRkV5ILplbd5Iwy798HPpxkN3AX8OvTKW3mDjUHhvJ/YjqKJPkwsAn46VnXMg1J3gJ8GviVGZcyC2sYLM18kMFvag8k+YmqemWWRU3BlcBnq+pPkrwf+FyS91TV92Zd2NHmSDpzH+W2Bd8fk2QNg1/bXpxKdZMz0u0akvws8Cngkqr6zpRqm7Rhc38b8B7g/iTPMFiL3NnBm6qjvOa7gZ1V9d2q+nfgawzC/mg2yry3ALcDVNW/ACcyuKlY78Z+25YjKdxHuW3BTuCq1r4MuLfauxFHsaHzTnIW8BcMgr2Htdc3HHTuVbWvqtZW1XxVzTN4v+GSqlqYTbljM8rP+t8zOGsnyVoGyzRPT7HGSRhl3v8BnA+Q5McZhPvSVKucjZ3AR9pVM+cC+6rq+VU946zfRT7AO8ZfY/CO+qda3x8y+AsNgxf6b4BF4MvAO2dd85Tm/U/AC8Aj7WvnrGue1tz3G3s/HVwtM+JrHgZLUk8CjwFXzLrmKc17I/AlBlfSPAL8/KxrHtO8Pw88D3yXwW9lW4CPAR9b9nrf2P5cHhvHz7m3H5CkDh1JyzKSpDEx3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KH/g/PtKk8R5qpAwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "draw_hist(yp[6])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD7CAYAAACRxdTpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAARrUlEQVR4nO3cf4xlZ33f8fcH20CaOLHBE2u7u2QdsoguUbN2pq6rVA0x+WGMmjVJaq2lBIPcbkAmP1RaYZKqIWmtmrbBKqprdYkdFpRgXJKIbTBpHGOEiGI7Y1ivvXYICyzybhbvBDA/hOLGzrd/3MflMp7Ze2fu3Lnjh/dLuppznvOce773zOxnz33Oj1QVkqS+PGfWBUiS1p/hLkkdMtwlqUOGuyR1yHCXpA4Z7pLUoZHhnuT5Se5L8kCSI0l+vbW/K8lnkxxqr92tPUnekeRoksNJLpryZ5AkLXHmGH2eAC6tqq8lOQv4WJIPtWX/tqrev6T/K4Gd7fWPgZvbT0nSBhkZ7jW4y+lrbfas9jrdnU97gHe39e5Jck6SLVV1cqUVzjvvvNqxY8f4VUuSuP/++/+6quaWWzbOkTtJzgDuB74PuKmq7k3yBuD6JP8euAu4rqqeALYCjw6tfry1rRjuO3bsYGFhYawPI0kaSPK5lZaNdUK1qp6qqt3ANuDiJN8PvAV4KfCPgBcAb15lUfuSLCRZWFxcXM2qkqQRVnW1TFU9DtwNXFZVJ2vgCeC3gYtbtxPA9qHVtrW2pe+1v6rmq2p+bm7ZbxWSpDUa52qZuSTntOlvA34M+IskW1pbgCuAh9oqB4HXtKtmLgG+fLrxdknS+htnzH0LcKCNuz8HuL2q/jDJh5PMAQEOAa9v/e8ALgeOAl8HXrfuVUuSTmucq2UOAxcu037pCv0LuHby0iRJa+UdqpLUIcNdkjpkuEtShwx3SerQWHeoSk/bcd0HZ7LdYze8aibblZ6tPHKXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdcibmJ6FZnUjkaRnD4/cJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0aGe5Jnp/kviQPJDmS5Ndb+wVJ7k1yNMn7kjy3tT+vzR9ty3dM+TNIkpYY58j9CeDSqvoBYDdwWZJLgLcBN1bV9wFfAq5p/a8BvtTab2z9JEkbaGS418DX2uxZ7VXApcD7W/sB4Io2vafN05a/IknWq2BJ0mhjjbknOSPJIeAUcCfwaeDxqnqydTkObG3TW4FHAdryLwMvXMeaJUkjjBXuVfVUVe0GtgEXAy+ddMNJ9iVZSLKwuLg46dtJkoas6mqZqnocuBv4J8A5SZ5+quQ24ESbPgFsB2jLvwv4wjLvtb+q5qtqfm5ubm3VS5KWNc7VMnNJzmnT3wb8GPAIg5D/mdbtauADbfpgm6ct/3BV1TrWLEkaYZznuW8BDiQ5g8F/BrdX1R8meRi4Lcl/BD4B3NL63wK8J8lR4IvA3inULUk6jZHhXlWHgQuXaf8Mg/H3pe1/A/yLdalOkrQm3qEqSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6NDLck2xPcneSh5McSfJLrf2tSU4kOdRelw+t85YkR5N8MslPTPMDSJKe6cwx+jwJvKmqPp7kbOD+JHe2ZTdW1X8d7pxkF7AXeBnw94E/SfKSqnpqPQuXJK1s5JF7VZ2sqo+36a8CjwBbT7PKHuC2qnqiqj4LHAUuXo9iJUnjWdWYe5IdwIXAva3pjUkOJ7k1ybmtbSvw6NBqxzn9fwaSpHU2drgn+Q7g94BfrqqvADcDLwZ2AyeB31zNhpPsS7KQZGFxcXE1q0qSRhgr3JOcxSDYf6eqfh+gqh6rqqeq6u+Ad/KNoZcTwPah1be1tm9SVfurar6q5ufm5ib5DJKkJca5WibALcAjVfX2ofYtQ91eDTzUpg8Ce5M8L8kFwE7gvvUrWZI0yjhXy/wQ8HPAg0kOtbZfAa5Kshso4Bjw8wBVdSTJ7cDDDK60udYrZSRpY40M96r6GJBlFt1xmnWuB66foC5J0gS8Q1WSOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHVoZLgn2Z7k7iQPJzmS5Jda+wuS3JnkU+3nua09Sd6R5GiSw0kumvaHkCR9s3GO3J8E3lRVu4BLgGuT7AKuA+6qqp3AXW0e4JXAzvbaB9y87lVLkk5rZLhX1cmq+nib/irwCLAV2AMcaN0OAFe06T3Au2vgHuCcJFvWu3BJ0spWNeaeZAdwIXAvcH5VnWyLPg+c36a3Ao8OrXa8tUmSNsjY4Z7kO4DfA365qr4yvKyqCqjVbDjJviQLSRYWFxdXs6okaYSxwj3JWQyC/Xeq6vdb82NPD7e0n6da+wlg+9Dq21rbN6mq/VU1X1Xzc3Nza61fkrSMca6WCXAL8EhVvX1o0UHg6jZ9NfCBofbXtKtmLgG+PDR8I0naAGeO0eeHgJ8DHkxyqLX9CnADcHuSa4DPAVe2ZXcAlwNHga8Dr1vPgiVJo40M96r6GJAVFr9imf4FXDthXZKkCXiHqiR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHRoZ7kluTnEry0FDbW5OcSHKovS4fWvaWJEeTfDLJT0yrcEnSysY5cn8XcNky7TdW1e72ugMgyS5gL/Cyts7/SHLGehUrSRrPyHCvqo8CXxzz/fYAt1XVE1X1WeAocPEE9UmS1mCSMfc3Jjnchm3ObW1bgUeH+hxvbZKkDbTWcL8ZeDGwGzgJ/OZq3yDJviQLSRYWFxfXWIYkaTlrCveqeqyqnqqqvwPeyTeGXk4A24e6bmtty73H/qqar6r5ubm5tZQhSVrBmsI9yZah2VcDT19JcxDYm+R5SS4AdgL3TVaiJGm1zhzVIcl7gZcD5yU5Dvwa8PIku4ECjgE/D1BVR5LcDjwMPAlcW1VPTaVySdKKRoZ7VV21TPMtp+l/PXD9JEVJkibjHaqS1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDI8M9ya1JTiV5aKjtBUnuTPKp9vPc1p4k70hyNMnhJBdNs3hJ0vLGOXJ/F3DZkrbrgLuqaidwV5sHeCWws732ATevT5mSpNU4c1SHqvpokh1LmvcAL2/TB4CPAG9u7e+uqgLuSXJOki1VdXLdKta3pB3XfXBm2z52w6tmtm1prdY65n7+UGB/Hji/TW8FHh3qd7y1SZI20MQnVNtReq12vST7kiwkWVhcXJy0DEnSkLWG+2NJtgC0n6da+wlg+1C/ba3tGapqf1XNV9X83NzcGsuQJC1nreF+ELi6TV8NfGCo/TXtqplLgC873i5JG2/kCdUk72Vw8vS8JMeBXwNuAG5Pcg3wOeDK1v0O4HLgKPB14HVTqFmSNMI4V8tctcKiVyzTt4BrJy1KkjQZ71CVpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdOnOSlZMcA74KPAU8WVXzSV4AvA/YARwDrqyqL01WpiRpNdbjyP1Hqmp3Vc23+euAu6pqJ3BXm5ckbaBpDMvsAQ606QPAFVPYhiTpNCYalgEK+OMkBfzPqtoPnF9VJ9vyzwPnT7iNTWvHdR+cdQmStKxJw/2fVtWJJN8N3JnkL4YXVlW14H+GJPuAfQAvetGLJixDkjRsomGZqjrRfp4C/gC4GHgsyRaA9vPUCuvur6r5qpqfm5ubpAxJ0hJrDvck357k7KengR8HHgIOAle3blcDH5i0SEnS6kwyLHM+8AdJnn6f362qP0ry58DtSa4BPgdcOXmZkqTVWHO4V9VngB9Ypv0LwCsmKUqSNBnvUJWkDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUoUkf+TtzPlNdkp7JI3dJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUoWf9pZCSNKlZXlJ97IZXTeV9PXKXpA4Z7pLUIYdlpBFm9ZV9Wl/X9a3BcJc2KR+toUlMbVgmyWVJPpnkaJLrprUdSdIzTSXck5wB3AS8EtgFXJVk1zS2JUl6pmkduV8MHK2qz1TV/wVuA/ZMaVuSpCWmFe5bgUeH5o+3NknSBpjZCdUk+4B9bfZrST45hc2cB/z1FN53UpuxLmsajzWNbzPWtelqytsmqul7VlowrXA/AWwfmt/W2v6/qtoP7J/S9gFIslBV89PcxlpsxrqsaTzWNL7NWNe3Uk3TGpb5c2BnkguSPBfYCxyc0rYkSUtM5ci9qp5M8kbg/wBnALdW1ZFpbEuS9ExTG3OvqjuAO6b1/mOa6rDPBDZjXdY0Hmsa32as61umplTVNN5XkjRDPjhMkjrURbiPetRBkn+d5OEkh5PclWTFy4c2sKbXJ3kwyaEkH9uIO3jHfSREkp9OUkk25KqCMfbVa5Mstn11KMm/nHVNrc+V7e/qSJLfnXVNSW4c2kd/meTxTVDTi5LcneQT7d/f5Zugpu9pOXA4yUeSbNuAmm5NcirJQyssT5J3tJoPJ7lo4o1W1bP6xeCE7aeB7wWeCzwA7FrS50eAv9em3wC8bxPU9J1D0z8J/NGsa2r9zgY+CtwDzG+S399rgf++yf6mdgKfAM5t898965qW9P8FBhcyzHo/7Qfe0KZ3Acc2QU3/C7i6TV8KvGcD/qb+GXAR8NAKyy8HPgQEuAS4d9Jt9nDkPvJRB1V1d1V9vc3ew+C6+1nX9JWh2W8Hpn3yY9xHQvwH4G3A30y5ntXWtZHGqelfATdV1ZcAqurUJqhp2FXAezdBTQV8Z5v+LuCvNkFNu4APt+m7l1m+7qrqo8AXT9NlD/DuGrgHOCfJlkm22UO4r/ZRB9cw+B9ymsaqKcm1ST4N/GfgF2ddU/squL2qNvJZs+P+/n66fV19f5Ltyyzf6JpeArwkyZ8muSfJZZugJmAw7ABcwDcCbJY1vRX42STHGVw99wuboKYHgJ9q068Gzk7ywinXNcq6P7Klh3AfW5KfBeaB/zLrWgCq6qaqejHwZuDfzbKWJM8B3g68aZZ1rOB/Azuq6h8CdwIHZlwPDC4j3gm8nMFR8juTnDPLgobsBd5fVU/NuhAG++ZdVbWNwdDDe9rf2iz9G+CHk3wC+GEGd89vhn21rma9k9fDyEcdACT5UeBXgZ+sqic2Q01DbgOumGZBjK7pbOD7gY8kOcZg3O/gBpxUHedRFV8Y+p39FvCDs66JwZHVwar626r6LPCXDMJ+ljU9bS/TH5KB8Wq6BrgdoKr+DHg+g+e7zKymqvqrqvqpqrqQQSZQVY9PsaZxrDYzRpv2iYQNOFFxJvAZBl9Dnz6B8rIlfS5kcJJl5yaqaefQ9D8HFmZd05L+H2FjTqiOs6+2DE2/GrhnE9R0GXCgTZ/H4Cv1C2f9+wNeChyj3cOyCfbTh4DXtul/wGDMfWq1jVnTecBz2vT1wG9Me1+1be1g5ROqr+KbT6jeN/H2NuJDbcBOu5zBkdOngV9tbb/B4Cgd4E+Ax4BD7XVwE9T034AjrZ67Txe0G1XTkr4bEu5j7qv/1PbVA21fvXQT1BQGw1gPAw8Ce2ddU5t/K3DDRvzextxPu4A/bb+7Q8CPb4Kafgb4VOvzW8DzNqCm9wIngb9l8K3vGuD1wOuH/p5uajU/uB7/9rxDVZI61MOYuyRpCcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QO/T89lDu7EB3RzgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "draw_hist(yp[7])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAANB0lEQVR4nO3df4zk9V3H8edLrtVEq1DveiGAHJpr9EyU4qWSaCxKUuFIerQmBBLtQdAzhvor1XjVP2hMGvEPbUKsmKslHI2FYrVypmglWEM0pXYRivyw7dkeclfKLdBiDUkVfPvHfknHY4+d2dnZ2X33+Ug2O/Od7873/clwz8x+d2ZIVSFJ6uVb5j2AJGntGXdJasi4S1JDxl2SGjLuktTQlnkPALB169basWPHvMeQpE3l/vvvf7qqti1324aI+44dO1hYWJj3GJK0qSR5/FS3eVpGkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGtoQ71CVpHnaceBjczv20Rsum8n9+sxdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkMrxj3JOUk+keTRJI8k+dVh+2uT3J3k88P3M4btSXJjkiNJHkpywawXIUn6/8Z55v4C8M6q2gVcCFyXZBdwALinqnYC9wzXAS4Fdg5f+4Gb1nxqSdIrWjHuVfVkVf3LcPlrwGPAWcBe4NCw2yHg8uHyXuDWWnIfcHqSM9d6cEnSqU10zj3JDuANwKeA7VX15HDTl4Htw+WzgCdGfuzYsO3k+9qfZCHJwuLi4qRzS5JewdhxT/IdwF8Av1ZV/zl6W1UVUJMcuKoOVtXuqtq9bdu2SX5UkrSCseKe5FUshf3Pquovh81PvXS6Zfh+Yth+HDhn5MfPHrZJktbJOK+WCfAB4LGq+sORmw4D+4bL+4A7R7a/fXjVzIXAcyOnbyRJ62DLGPv8GPBzwL8meXDY9tvADcAdSa4FHgeuGG67C9gDHAGeB65Zy4ElSStbMe5V9Y9ATnHzxcvsX8B1U84lSZqC71CVpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkMrxj3JzUlOJHl4ZNu7kxxP8uDwtWfktnclOZLks0l+elaDS5JObZxn7rcAlyyz/b1Vdf7wdRdAkl3AlcAPDj/zx0lOW6thJUnjWTHuVXUv8OyY97cXuL2qvl5VXwSOAG+cYj5J0ipMc879HUkeGk7bnDFsOwt4YmSfY8O2l0myP8lCkoXFxcUpxpAknWy1cb8J+D7gfOBJ4A8mvYOqOlhVu6tq97Zt21Y5hiRpOauKe1U9VVUvVtX/Au/nG6dejgPnjOx69rBNkrSOVhX3JGeOXH0r8NIraQ4DVyb51iTnATuBf55uREnSpLastEOS24CLgK1JjgHXAxclOR8o4CjwiwBV9UiSO4BHgReA66rqxZlMLkk6pRXjXlVXLbP5A6+w/3uA90wzlCRpOr5DVZIaMu6S1JBxl6SGjLskNbTiH1R1ajsOfGwuxz16w2VzOa6kzcNn7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGtr0b2Ka1xuJJGkj85m7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktTQinFPcnOSE0keHtn22iR3J/n88P2MYXuS3JjkSJKHklwwy+ElScsb55n7LcAlJ207ANxTVTuBe4brAJcCO4ev/cBNazOmJGkSK8a9qu4Fnj1p817g0HD5EHD5yPZba8l9wOlJzlyjWSVJY1rtOfftVfXkcPnLwPbh8lnAEyP7HRu2SZLW0dR/UK2qAmrSn0uyP8lCkoXFxcVpx5AkjVht3J966XTL8P3EsP04cM7IfmcP216mqg5W1e6q2r1t27ZVjiFJWs5q434Y2Ddc3gfcObL97cOrZi4Enhs5fSNJWidbVtohyW3ARcDWJMeA64EbgDuSXAs8Dlwx7H4XsAc4AjwPXDODmSVJK1gx7lV11SluuniZfQu4btqhJEnT8R2qktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpoS3zHkDS8nYc+Nhcjnv0hsvmclyY35o7Mu7SCgyONiNPy0hSQ8Zdkhoy7pLUkHGXpIaMuyQ15KtlNJFvxpfnfbPx1UE9+Mxdkhoy7pLU0FSnZZIcBb4GvAi8UFW7k7wW+DCwAzgKXFFVX5luTEnSJNbinPtPVtXTI9cPAPdU1Q1JDgzXf2sNjqOB50QlrWQWp2X2AoeGy4eAy2dwDEnSK5g27gX8XZL7k+wftm2vqieHy18Gti/3g0n2J1lIsrC4uDjlGJKkUdOelvnxqjqe5HXA3Un+bfTGqqoktdwPVtVB4CDA7t27l91HkrQ6Uz1zr6rjw/cTwEeBNwJPJTkTYPh+YtohJUmTWXXck3x7kte8dBl4M/AwcBjYN+y2D7hz2iElSZOZ5rTMduCjSV66nw9V1d8m+TRwR5JrgceBK6YfU5I0iVXHvaq+APzwMtufAS6eZihJ0nR8h6okNWTcJakh4y5JDRl3SWrIuEtSQ/7POrQp+GFp0mR85i5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqaGZxT3JJUk+m+RIkgOzOo4k6eVmEvckpwHvAy4FdgFXJdk1i2NJkl5uVs/c3wgcqaovVNV/A7cDe2d0LEnSSbbM6H7PAp4YuX4M+NHRHZLsB/YPV/8ryWdnNMt62go8Pe8hZsB1bS6uaxPJ70+1rnNPdcOs4r6iqjoIHJzX8WchyUJV7Z73HGvNdW0urmtzmdW6ZnVa5jhwzsj1s4dtkqR1MKu4fxrYmeS8JK8GrgQOz+hYkqSTzOS0TFW9kOQdwMeB04Cbq+qRWRxrg2l1mmmE69pcXNfmMpN1papmcb+SpDnyHaqS1JBxl6SGjPsqrPTRCkmuTrKY5MHh6+fnMeekxvnIiCRXJHk0ySNJPrTeM67GGI/Xe0ceq88l+eocxpzYGOv6niSfSPJAkoeS7JnHnJMYY03nJrlnWM8/JDl7HnNOKsnNSU4kefgUtyfJjcO6H0pywdQHrSq/Jvhi6Q/E/w58L/Bq4DPArpP2uRr4o3nPOoN17QQeAM4Yrr9u3nOvxbpO2v+XWXoBwNxnX4PH6yDwS8PlXcDRec+9Bmv6c2DfcPmngA/Oe+4x1/YTwAXAw6e4fQ/wN0CAC4FPTXtMn7lPrutHK4yzrl8A3ldVXwGoqhPrPONqTPp4XQXcti6TTWecdRXwncPl7wK+tI7zrcY4a9oF/P1w+RPL3L4hVdW9wLOvsMte4NZach9wepIzpzmmcZ/cch+tcNYy+/3M8OvVR5Kcs8ztG80463o98Pok/5TkviSXrNt0qzfu40WSc4Hz+EY8NrJx1vVu4GeTHAPuYum3ko1snDV9BnjbcPmtwGuSfPc6zDZrY/93Oi7jPht/Deyoqh8C7gYOzXmetbKFpVMzF7H0DPf9SU6f50Br7ErgI1X14rwHWSNXAbdU1dks/dr/wSSb/d/8bwBvSvIA8CaW3vne5fFaU5v9gZ6HFT9aoaqeqaqvD1f/FPiRdZptGuN8ZMQx4HBV/U9VfRH4HEux38gm+SiMK9kcp2RgvHVdC9wBUFWfBL6NpQ/f2qjG+bf1pap6W1W9AfidYdtX123C2Vnzj2wx7pNb8aMVTjpX9hbgsXWcb7XG+ciIv2LpWTtJtrJ0muYL6zjjaoz1URhJvh84A/jkOs+3WuOs6z+AiwGS/ABLcV9c1yknM86/ra0jv328C7h5nWeclcPA24dXzVwIPFdVT05zh3P7VMjNqk7x0QpJfhdYqKrDwK8keQvwAkt/RLl6bgOPacx1fRx4c5JHWfpV+Der6pn5Tb2yMdcFSyG5vYaXLmx0Y67rnSydOvt1lv64evVGXt+Ya7oI+L0kBdwLXDe3gSeQ5DaWZt86/A3keuBVAFX1Jyz9TWQPcAR4Hrhm6mNu4MdakrRKnpaRpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGvo/zwMX4X4CHFUAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "draw_hist(yp[8])"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "pytorch",
   "language": "python",
   "name": "pytorch"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
